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ABSTRACT 

On-going projects to discover Type la supernovae at redshifts z ~ 0.3 — 1, 
coupled with improved techniques to narrow the dispersion in SN la peak mag- 
nitudes, have renewed the prospects for determining the cosmic deceleration pa- 
rameter q . We estimate the expected uncertainty in the Hubble diagram deter- 
mination of go due to weak lensing by structure in the universe, which stochas- 
tically shifts the apparent brightness of distant standard candles. Although the 
results are sensitive to the density power spectrum on small scales, the induced 
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flux dispersion a m < 0.04f2„; mag for sources at z < 0.5, well below the "in- 
trinsic" spread of nearby SN la magnitudes, om — 0.2 mag. Thus, density 
inhomogeneities do not significantly impact the current program to measure go, 
in contrast to a recent claim. If, however, light-curve shape and other calibrators 
can reduce the effective intrinsic spread to 0.1 mag at high z, then weak lensing 
could increase the observed spread by 30 % in an Q m ,o = 1 universe for SNe at 



Subject headings: cosmology: large-scale structure of the universe, supernovae 

1. Introduction 

The determination of the cosmological parameters via the "classical" cosmological tests, 
such as the redshift-magnitude relation (Hubble diagram), remains a holy grail for observa- 
tional astronomy. Past attempts to measure the deceleration parameter go using galaxies as 
standard candles have foundered on the uncertainty in galaxy luminosity evolution ( [lmsley 
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1972, [Peebles 1993Q . Recently, there has been renewed hope that g can be determined, 
thereby constraining the mean mass density and the cosmological constant, by using distant 
Type la supernovae (e.g., |Goobar fc Perlmutter 1995]). 



Progress has come on two fronts. First, there is growing evidence that samples of SNe 
la, when suitably culled to excise peculiar lightcurves or spectra and cases of significant 
host-galaxy extinction, may provide reasonably good standard candles, with a dispersion 
om ~ 0.3 in peak absolute magnitude. Moreover, using the observed correlation between 
lightcurve shape and peak luminosity ( |Hamuy etal. 1995] , |Riess, Press, fc Kirshner 1995| ) as 



well as spectroscopic features that have been found to correlate with luminosity ( [Branch 



[Nugent, fc Fisher 1996| , [Nugent, etal. 1995] ), the effective dispersion of these 'standardized 



candles' can apparently be reduced to cr M ~ 0.1 — 0.2 mag. Second, several groups have 
begun observing campaigns to discover a large number of high-redshift la supernovae, with 
coordinated follow-up programs to measure the lightcurves of SN candidates; more than 20 
candidates at z ~ 0.4 — 0.6 have been found so far (e.g., |Perlmutter, etal. 1995| , [Schmidt 
etal. 1995| , [Perlmutter, etal. 1996j ). 



The case for SNe la as a population of 'standard' candles is also on reasonable physical 
footing. While the details of the explosion mechanism remain poorly understood, there is 
general consensus that SNe la are thermonuclear explosions of accreting Carbon-Oxygen 
white dwarfs. On the other hand, it is not yet clear whether the white dwarf progenitors 
must be near the Chandrasekhar mass (e.g., |Hoeflich etal. 1996|) or not ( [Livne fc Arnett 



1995|) , and thus the theoretical intrinsic spread in SN la luminosity is still a matter of some 
debate. In either case, there is optimism that, based on the correlations above, SNe la 
can be used to measure extragalactic distances and thereby to determine the cosmological 
parameters. 

Recently, however, [Kantowski, Vaughan, fc Branch (1995)| have argued that large-scale 



structure provides a stumbling block to the accurate determination of go from standard can- 
dles. A bundle of light rays from a distant source is sheared and focused due to deflection by 
intervening density inhomogeneities. Since gravitational lensing conserves surface brightness, 
the change in the cross-section of the bundle implies that the image of the source appears 
magnified (or demagnified) relative to a homogeneous universe. Using the "swiss-cheese" 
model of large-scale structure, Kantowski etal. find that the resulting change in received 
flux (apparent magnitude) could lead to a systematic underestimate of g if one interprets 
the observations assuming a homogeneous universe. For example, for a source at z = 0.458 
(the redshift of SN 1992bi ( [Perlmutter, etal. 1995[ )), they find that the resulting error in go 
can be as large as Sqo ~ — 0.33go. This claim is consonant with other studies which found 
that large-scale structure can significantly change the proportionality between angular-size 
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distance and redshift, and therefore lead to difficulties, e.g., for the determination of H from 
gravitational lens time delays ( |Kantowski 1969] , Pyer fc Roeder 1972| , |Alcock fc Anderson 



19851 , 1986, |Watanabe etal. 1992^ ^asaki 1993| ). 



In this Comment, we reevaluate this issue by estimating the rms fluctuation in the am- 
plification of distant sources in a perturbed Friedmann-Robertson- Walker (FRW) universe. 
We find that the expected effects are smaller than those found by Kantowski, etal., of order 
several percent at most for sources at z < 0.5, and subdominant in comparison to the ~ 20 
% intrinsic spread of nearby SN la magnitudes. Moreover, flux conservation implies that 
the magnification shift is random, with zero mean over all lines of sight flWeinberg 1976| ), 
not a systematic offset. As a result, the amplification due to large-scale structure will not 
seriously impact the accuracy of go measurements in current surveys. 

The primary reason for this different conclusion is that the swiss-cheese model does 
not conform to our current understanding of the large-scale mass distribution of the uni- 
verse. In particular, it does not accurately reflect the observational information gained from 
galaxy redshift and peculiar velocity surveys and the cosmic microwave background (CMB) 
anisotropy. This point has been made recently by |Seljak (1994)| and Bar-kana (1995) in the 
context of the determination of Ho from QSO lens time delays. The argument that these 
effects should be small has also been made quantitatively by Peebles (1993). The issue was 
first laid out clearly by |Gunn (1967)| , who showed that the rms fluctuation in the apparent 
brightness of a distant source can be expressed as a radial integral of the two-point correlation 
function of the mass distribution. Gunn's argument was updated to reflect the substantial 
advances in our understanding of large-scale structure in more recent numerical ( |Jaroszynski, 



etal. 199"0p and analytic (Pabul fc Lee 199~1|) studies (unbeknownst to the present author 



until this work was completed). In fact, Babul Sz Lee's 'DP' model is close to the model for 
the power spectrum that we adopt below, although there are some important quantitative 
differences which we outline later. 

From the smallness of the CMB anisotropy on large scales to observations of galaxies 
and galaxy clusters on small scales, we have strong indications that the spacetime metric of 
the universe is well- described by a weakly perturbed FRW model. In the longitudinal gauge, 
the line element can be written 



ds 2 



a 2 (r) 



[1 + 2<p)dr 2 + (1 - 2<p)[d X 2 + F 2 ( X )(d6 2 + sin 2 6d<p 2 )} 



where a(r) is the cosmic scale factor, \ denotes the comoving radial coordinate, and r = 
J dt/a is the conformal time, with tq denoting the present. The function F(x) depends 
on the spatial curvature K: F(x) = K~ 1 ' 2 sin K l ' 2 x for K > 0, F = x for K = 0, and 
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F = (— ii') -1 / 2 sinh(— Ky^x for X < 0. The curvature can be expressed in terms of the 
present density parameter and the Hubble parameter H = H(t ), K = (Q — 1)HqQq, 
where Q = Q m + Q\ includes both non-relativistic matter (m) and vacuum energy density 
(the cosmological constant A) and O m = p/p c = 8irGp/3H 2 , with p(r) the mean density 
of matter. The metric perturbation variable <fi is the relativistic analog of the Newtonian 
gravitational potential; over scales less than the Hubble length H^ 1 it obeys the Poisson 
equation, 



V 2 = hl m H 2 a 2 5 (2) 

where the density contrast <5(x, r) = p(x, r)/p — 1. Observationally, the fluctuations in 
spatial geometry correspond to < 1CT 4 — 1CT 5 from galaxy scales to the Hubble radius. 
We will assume 1 but do not place any restrictions on the amplitude of 6. 

For a photon with direction n, the null geodesic equation in the metric (JI|) determines 
the rate of change in propagation direction due to inhomogeneities along the light path; to 
lowest order in 0, the geodesic can be parametrized by the radial coordinate Xi yielding 
dh/dx = — 2Vj_0, where the derivative is transverse to the line of sight; we are implicitly 
assuming small deflection angles, so that the two-dimensional sphere can be approximated 
by a plane perpendicular to the unperturbed line of sight. Directions in this plane can be 
parameterized by the two-dimensional angle 0. A photon observed at direction 6 would 
have been seen at direction f3 = + 66 in the absence of weak lensing effects, i.e., 86 is 
the net transverse deflection of the light rays due to density fluctuations. From the geodesic 
equation, for a source at comoving distance Xs we have (e.g., |Kaiser 1992| , pr*yne fc Birkinshaw 
1996| , |Villumsen 1995| , |Seljak 1996|) , defining x — ~ n X> 



66 = -2 £ 3 g£frX> VL0(x, r = r - X )dx ■ (3) 

Here we have implicitly replaced the perturbed by the unperturbed path in the integral; 
while this is not a good approximation (e.g., Frieman, Harari, fc Surpi 1994 ), we will only 



make use of the weaker assumption that the statistical properties of the potential field are 
identical along the perturbed and unperturbed paths (e.g., [Kaiser 1992|) . 

The amplification of the observed image relative to the (unlensed) source is given by 
A = 1/det M, where the 2x2 amplification matrix iWy = d(3i/d9j = 5ij + and 
$jj = d59i/d9j. We can decompose the deformation tensor as $y = — k8^ + 7^, where the 
trace (k), the expansion, describes the uniform dilation or contraction of ray bundles, and the 
traceless part (7^) describes their shear. In the limit of small deflection angles, we thus have 



- 5 - 



A = 1/[(1 — ft) 2 — 7 2 ] — 1 + 2k, where ±7 are the eigenvalues of 7^ and k = — ($n + $22)/2. 
The flux perturbation is therefore 8 A = A — 1 = — Tr$. Using eqns. (fj, |3|), we find 

5 A = 3(g °;; )2 f-° f ^(X)i^ - r = r - X ) (4) 

^(Xs) Jo o(r - x) 

where f2 m) o is the present matter density parameter. 

To find the rms flux amplification along a random line of sight, we assume <5(x, r) can 
be described as a continuous, homogeneous random process; this ignores discreteness in the 
mass density, an excellent approximation since the dark matter very likely consists of objects 
of mass less than 10 6 M Q . It is convenient to Fourier-transform the density field, S(x) = 
(27r) -3 / d 3 k 5(k) exp(ik • where the flat-space transform suffices for the small angles we 
are considering. Defining the density power spectrum via (5(k)<5*(k')) = (27r) 3 P(k)<5|)(k— k'), 
we have (defining a\ = ((5 A) 2 )) 

a> A = ffcGUW f ^)^*>,l + ,(r)T- J (5) 

where A 2 (A;, z) = k 3 P(k, z)/2n 2 = da 2 /dink is the contribution per logarithmic wavenumber 
interval to the variance of the density field. Here, z(t) is the redshift at epoch r = r — 
given implicitly by 

, , 1 r z dz' 

x(z) = / (6) 

H a Jo [n m>0 (l + z'f + (1 - Q mfi - Q A )(1 + z') 2 + Q A ] 1/2 

For Q A = 0, this has a well-known analytic solution for z(x) ( [Mattig 1958| ). In (|^), we have 
used a small-angle approximation, so that only waves nearly perpendicular to the line of 
sight contribute to the amplification ( Blandford, etal. 1991| , [Kaiser 1992j ); in addition, we 



have assumed that the density correlation length is small compared to the Hubble radius, so 
we only include contributions to (5(k, r)<5*(k', r')) from equal times, r = t'. 

Over a limited range of wavenumber k and redshift z, the power spectrum scales with 
time as A 2 (k,z) = A 2 (k)(l + z)~ e , separating the integrals in (|^). At small fc, A < 1, 
and it obeys linear perturbation theory, e = 2. At large k, where A ^> 1, the clustering is 
highly non-linear and has approximately reached virial equilibrium, i.e., galaxies and groups 
of galaxies reach a fixed physical size; in this 'stable clustering' regime, if the two-point 
density correlation function obeys £(r) ~ r~ 7 , then e = 3 — 7. Since the observed galaxy 
correlation function is a power-law with 7 ~ 1.8 at r < lOh^ 1 Mpc, we expect e ~ 1.2 on 
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small scales, although with the caveat that the galaxy and mass distributions may differ 
significantly on these scales (see below). A third more radical possibility is that the observed 
clustering is merely "painted on" and expands with the Hubble flow; in this case, e = 0. 
For the numerical estimates below, we will thus consider the range e = — 2. Since the rms 
amplification is dominated by structure in the non-linear regime, we expect e ~ 1.2 to be 
the most accurate representation of the evolution. Recent N-body simulations confirm that 
the growth factor is intermediate between the stable clustering and linear regimes on small 
scales at recent epochs flColin, Carlberg, fc Couchman 1996Q , while it may be faster than 
linear (e > 2) on intermediate scales; in either case, the results presented below constitute 
an upper limit on the lensing effect. 

To model the present power-spectrum on small scales, we could simply Fourier trans- 
form the galaxy correlation function, £ 9 (r) = (r/ro)~ L8 , where the galaxy correlation length 
ro i9 — 5Ah~ l Mpc. However, this does not take into account the bias between the galaxy 
and mass distributions, i.e., the likelihood that light does not trace mass on these scales. 
To remedy this, we can incorporate dynamical information. From the (simplified version 
of the) cosmic virial theorem (Peebles 1980, 1993), the predicted pairwise velocity disper- 
sion on small scales is roughly cr^(r) ~ [3/2(3 — 7)]-£/~Qfi mi0 r 2 (r /r) 7 . Assuming the matter 
correlation function has the same slope as that of the galaxies, this yields an estimate for 
the density correlation length, r ~ 5.4/i -1 Mpc(0.1/fi mi0 )°- 55 (cr t ,(l/i _1 Mpc)/300 km/sec) 11 . 
The standard estimates of the galaxy pairwise velocity dispersion at r ~ \hr x Mpc sep- 
aration have been a v ~ 300 km/sec ( Pavis fc Peebles 1983| , |Bean, etal. 19831) . Recent 



redshift surveys, however, have yielded higher values, a v ~ 500 — 700 km/sec flGuzzo, etal 



1995| , Marzke, etal. 1995| ). Moreover, Somerville, Davis, & Primack (1996) and Somerville, 



Primack, & Nolthenius (1996) have found that estimates of cr„(r) in both simulations and 
galaxy catalogs show large scatter. However, when the cores of rich clusters, which contain 
only a small fraction of the mass, are excluded, it appears that a v ~ 300 km/sec. In any 
case, to obtain an upper bound for the weak lensing effect, we will take o v {\hr x Mpc) = 650 
km/sec, implying A 2 (k) ~ 8.7(k/h Mpc _1 ) L8 ^~* (cr,;(l)/650) 2 . Note that the assumption 
made here that the galaxy and density correlation functions have the same small-scale slope 
is not well motivated on highly non-linear scales, and the results below should be interpreted 
with this caveat. 

For this model of the power spectrum, and in general for 7 > 1, the /c-integral in (^) 
diverges in the ultraviolet and must be cut-off. Physically, the slope of the density correlation 
function must fall below unity below some length scale (even though the galaxy correlation 
function is an unbroken 7 ~ 1.8 power law down to r ~ lOh^ 1 kpc). This flattening is 
expected to happen at least at the scale of individual galaxy halos ( |Kaiser 1992j ): below 
this scale, £(r) is dominated by correlations within individual halos, yielding £(r) ~ r -(2f-3) 
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for halos with density profile p(r) ~ r~ v ( [Peebles 1971 |Mc(Jlelland fc Silk 197^ , |Sheth fc 



Jain 1996| ); for nearly isothermal halos, v ~ 2, the corresponding slope is 7 ~ 1. Up to 



logarithmic corrections, we model this effect by imposing a cutoff in ([5]) at the approximate 
halo scale, k c = 1 / kpc). Using this model in (^), the rms flux perturbation at fixed 
z s scales approximately as flH 2 (a v (1)/650)(/c c /10/i Mpc -1 ) 0,4 . 

Fig.l shows the dispersion in apparent magnitude for distant standard candles as a 
function of redshift, a m = 1.086a a mag, for cosmological models with Q m ,o = 0.1, 0.3, and 
1. For fl m fl = 1, we show results for e = 0, 1.2, and 2 to bracket the plausible range of 
evolution models. For the other cases, we show results for e = 1.2 only; since structure 
formation freezes out early in f2 m < 1 models (at 1 + Zf — ^^o)> stable clustering should 
be a reliable prescription here. For Q m ,o — 0.3, we show results for a spatially flat model 
with non-zero cosmological constant, Q\ = 0.7, and for an open model with A = 0. 

Since the imposition of a sharp cutoff in the density power spectrum appears to be a 
rather crude approximation, as a check we have also explored models for the small-scale 
power spectrum based on the cold dark matter (CDM) model and CDM with non-zero 
A, with scale-invariant primordial fluctuations, extended into the non- linear regime using 
scaling formulae derived from N-body simulations ( Hamilton etal. 1991 , Peacock fc Dodds 



1991 Uain, Mo, fc White 1995| , [Baugh fc Gaztanaga 1996| , Peacock fc Dodds 1996|) . In these 



models, P{k) scales as k at small wavenumber but turns over to k~ 3 at large k, yielding 
much better convergence for the flux dispersion a a- We normalize such models by requiring 
that they reproduce the observed galaxy cluster X-ray temperature distribution function 
according to the predictions of Press-Schechter theory ([White, Efstathiou, fc Frenk 1993| ); 
for CDM models, this corresponds approximately to imposing the constraint that the linear 
theory rms mass fluctuation in a sphere of radius 8 h^ 1 Mpc is <Tg = OfiQ^Q^ 1 " 1 ' , where 
C = 0.36 + 0.3in rnfl - 0.28^ for open models (A = 0) and C = 0.59 - 0.16fi mj0 + 0.06fl^ 
for spatially flat (non-zero A) models ( |Viana fc Liddle 1995j ). The results for the amplification 
dispersion in these cluster-normalized CDM models are shown in Fig. 2 for the same model 
parameters as in Fig. 1. The results for a a agree reasonably well with those of the power- law 
model, although they are somewhat higher for the models with Q m fi < 1- This difference 
is due in part to the fact that the cluster normalization yields more small-scale power in 
these models than the cosmic virial theorem normalization for the corresponding power-law 
model. Thus, while the uncertainties in a v and in the accuracy of the cosmic virial theorem 
are still significant, this comparison suggests that our estimate for a a should be accurate to 
within a factor of two. 

The implications of weak lensing for the determination of go follow directly from Figs. 
1 and 2. For sources at z < 0.7, a m < 0.06, well below the "intrinsic" 0.2 — 0.3 mag spread 
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in nearby SN la magnitudes. At z <C 1, from the redshift-magnitude relation, the resulting 
'lcr' uncertainty in the deceleration parameter, q = fi mj0 /2 — Q^, for a single source at 

1 /2 

redshift 2; is o" go ~ (Ta/z. For example, for z = 0.5 and A = 0, we find a go ~ 0.1g 
(< 0.07 for flmfi < 1). In the future, if SN la searches discover sources at z > 1, then weak 
lensing may be a significant factor. In particular, if light-curve shape and other correlators 
with peak magnitude can reduce the effective intrinsic spread to 0.1 mag (even at high z), 
then density fluctuations could increase the observed dispersion in an fi m0 = 1 universe for 
sources at z = 1 by of order 30%. Since the amplification is caused by the foreground mass 
distribution, one could in principle use the angular correlation function of 8 A to probe the 
large-scale mass power spectrum; in practice, this will require thousands of well-measured 
SNe la at redshifts z > 0.5 spread over hundreds of square degrees. This may be possible 
with the Next Generation Space Telescope. 

The results shown here can be compared to those of Babul & Lee (1991). While they 
considered only the Einstein-de Sitter fl m ,o = 1 case, we have presented results for arbitrary 
fl m and A. Quantitatively, their 'DP' model for the power spectrum is quite similar to 
that used here, but it overestimated the rms amplification a a because it did not include the 
constraint from galaxy pairwise velocities on small scales. They also presented results for the 
Q = 1 CDM model with bias factor b ~ 2.5; this normalization, which reflected earlier, lower 
estimates of the small-scale pairwise velocity dispersion, is now known to be unacceptably 
low and therefore substantially underestimates a a- Thus, while their results bracket those 
here, the estimates in Figs. 1 and 2 should be more accurate. 

A final issue of concern for these surveys is amplification bias. A fraction of the SNe in a 
magnitude-limited survey are strongly amplified (5 A ^> a a) by foreground mass concentra- 
tions very close to the line of sight; this would cause large systematic errors in the estimate of 
go- However, the lensing galaxy or cluster responsible could generally be seen and such events 
removed from the sample. Even if the lenses are too faint to be detected, the probability for 
such strong lensing events at moderate redshift, z s < 1, is known to be very small, based 
on the low incidence of multiply imaged QSOs. While the optical depth r for significant 
amplification by foreground galaxies (e.g., 5 A > 0.1) is much higher than that for multiple 
imaging, it is still quite low. For example, for Q m ,o — 1 and assuming isothermal galaxy 
halos extend to r > 50/i _1 kpc, integration over the galaxy luminosity function with the 
Faber- Jackson relation gives t(5A > 0.1, z s — 1) ~ 6.5 x 10~ 3 for amplification by individual 
foreground galaxies; for a non-zero cosmological constant satisfying = 1 — Q m < 0.8, 
r(> 0.1, 1) < 0.02. The result of these near encounters is a very small non-Gaussian tail 
in the amplification distribution at large SA. On the other hand, if a substantial fraction 
of the cosmic density is in compact objects, the high amplification tail can be much more 
significant ( [Schneider fc Wagoner T9871 , ILinder, Schneider, fc Wagoner 19881 , |R aU ch 1991| ). If 
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high amplification (> 1 mag) events are not soon found in the high-redshift SN la searches, 
one will be able to constrain the contribution to Q from objects of mass > 1O _3 M . 

For completeness, I note that weak lensing does not affect the shape of SN la lightcurves: 
for a source at redshift z s , the fall-off from the peak is given by the usual time dilation factor, 
At = Aij(l + z s ), where At, is the fall-off timescale in the SN rest frame. 

Note added: After this work was accepted for publication, two recent numerical efforts 
have clarified the weak lensing effects of large-scale structure ( |Wambsganss, etal. 1996 , 
D. E. Holz and R. M. Wald, to be published). Using ray-shooting techniques through N- 
body and phenomenologically constructed universes, it has been directly confirmed that the 
mean flux from a distant standard candle is that given by the standard FRW luminosity 
distance. However, the distribution of the amplification is not symmetric about this mean. 
The majority of lines of sight pass through regions underdense compared to the mean, 
leading to modest de-amplification, while a smaller number of rays pass near dense mass 
concentrations and suffer substantial amplification. Thus, for a survey with a finite number 
of sources, there can be a small bias in the results for q , but the effect is very small at the 
current redshifts (z < 0.6) probed by the SN la surveys. 

I thank D. Holz, A. Olinto, S. Perlmutter, U. Seljak, J. Silk, A. Stebbins, and R. Wald 
for conversations, and the Institute for Nuclear and Particle Astrophysics (LBL) and the 
Center for Particle Astrophysics (UC Berkeley) for hospitality while this work was being 
completed. This research was supported in part by the DOE and by NASA grant NAG5- 
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Fig. 1. — Fig. 1. Dispersion in flux (in magnitudes) for a source at redshift z, for fi m)0 = 0.1, 
0.3, 1.0. Solid curves assume stable clustering, e = 1.2. For f2 m o = 1, dashed curve 
corresponds to linear theory evolution, e = 2, and dotted curve corresponds to 'painted on' 
structure (e = 0). Lower curves marked fi mj o = 0.1,0.3 are open models; curve marked 
'A' is a flat universe with Q\ = 0.7. The small-scale power spectrum has been cut off at 
k c = l/(100/i _1 kpc); the dispersion scales as a oc k® A . 
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Fig. 2. — Fig. 2. Flux dispersion vs. redshift for models with the same cosmological 
parameters as in Fig. 1, but for the CDM model of structure formation instead of the 
phenomenological power-law model. The models are standard CDM n m0 — 1, h — 0.5 
(solid), a flat ACDM model with JIa = 0.7, h — 0.7 (dotted), and an open CDM model with 
^m,o = 0.3, h = 0.7 (dashed). 



